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ABSTRACT 

The magnification of distant sources by mass clumps at lower (z < 1) red- 
shifts is calculated analytically. The clumps are initially assumed to be galaxy 
group isothermal spheres with properties inferred from an extensive survey. The 
average effect, which includes strong lensing, is exactly counteracted by the beam 
divergence in between clumps (more precisely, the average reciprocal magnifica- 
tion cancels the inverse Dyer-Roeder demagnification). This conclusion is inde- 
pendent of the matter density function within each clump, and remains valid 
for arbitrary values of Q m and Q\. When tested against the cosmic microwave 
background data, a rather large lensing induced dispersion in the angular size 
of the primary acoustic peaks of the TT power spectrum is inconsistent with 
WMAP observations. The situation is unchanged by the use of NFW profiles 
for the density distribution of groups, which led in fact to slightly larger fluctu- 
ations. Finally, our formulae are applied to an ensemble of NFW mass clumps 
or isothermal spheres having the parameters of galaxy clusters. The acoustic 
peak size dispersion remains unobservably large, and is also excluded by WMAP. 
For galaxy groups, two possible ways of reconciling with the data are proposed, 
both exploiting maximally the uncertainties in our knowledge of group proper- 
ties. The same escape routes are not available in the case of clusters, however, 
because their properties are well understood. Here we have a more robust con- 
clusion: neither the NFW profile nor isothermal sphere profiles are an accurate 
description of clusters, or important elements of physics responsible for shaping 
zero curvature space are missing from the standard cosmological model. When 
all the effects are accrued, it is difficult to understand how WMAP could reveal 
no evidence whatsoever of lensing by groups and clusters. 

1. Introduction - observed statistics of galaxy groups 
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Light as it propagates through the near Universe will encounter inhomogeneities. The 
motivation which initiated the present paper is to investigate the role played by groups of 
galaxies in the phenomenon of lensing and global curvature, by using (perhaps for the first 
time) real observational data. Recently, large databases on groups have become available, 
including in particular an ESO survey (Ramella et al 2002), from which some general prop- 
erties of 1,168 groups can be inferred. Specifically concerning the lensing performance of the 
groups, one needs to know their mass and velocity dispersion. These are shown in Figure 1. 
We then found the mean virial mass per group to be 

M group w 1.15 x 1O 14 M . (1) 

The datum can then be used to estimate the value of fig roU p, since these 1,168 systems were 
identified within a surveyed volume of 0.0075 Gpc 3 (comoving volume of a pie subtending 
4.69 sr at the observer and of radius extending to redshift z = 0.04, in a h = 0.71, VL m = 
0.27, Ha = 0.73 cosmology). Hence the number density of groups is: 

n gr0 up = n = 1.56 x 10~ 4 Mpc -3 , (2) 

and is in fact a lower limit, because no attempt was made here to correct n groU p for selection 
effects. Eqs. (1) and (2) point to a mass density of 

Mgroup ~ 0.135. (3) 

Given that most galaxies exist in groups, the value is in agreement with the conclusion 
of Fukugita (2003) and Fukugita, Hogan, & Peebles (1998), who found (after careful mass 
budget accountancy) that at low redshifts matter amounting to 50 % the total density of 
Q m ~ 0.27 is connected with galaxies and galactic environments. 

How does such a form of mass concentrations affect the propagation of light? To answer 
this question, some information about the matter profile in groups is necessary. The limited 
isothermal sphere model, wherein the internal matter density falls radially as 1/r 2 to some 
cutoff radius R, is sometimes advocated (Zabludoff & Mulchaey 1998, Mulchaey 2000). In 
this model the value of R is related to the total mass M of the group by the equation 

~R~ = 2<T ' (4) 

where a is the dispersion velocity of the group. Since the median dispersion velocity of the 
ESO group sample is a ~ 270 km s _1 (Ramella et al 2002), Eq. (1) and (3) may be coupled 
to provide an estimate of the cutoff radius R as R « 3 Mpc. 



2. Weak lensing: direct incorporation of galaxy group data 
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Concerning the observed breakdown of the Universe's total mass density, viz. Q,\ = 
0.73 dark energy and Q m = 0.27 matter (Bennett et al 2003), while at recent epochs (z ~ 1 
or less) the A component may remain smooth, the matter is certainly known to be clumped 
into mass concentrations, with galaxy groups forming an important subclass. 

The efforts to date on the global (all sky) influence exerted by weak lensing involve 
primarily N-body simulations (e.g. Wambsganss et al 1997, Barber 2000), complemented 
by some development on theoretical methods (Dalai et al 2003). There are, however, two 
areas of neglect: (a) if groups were adequately represented in previous works, it is doubtful 
whether direct observational properties were involved; (b) analytical formulae are generally 
lacking, even though they should be provided to the furthest extent possible - the physics is 
always much clearer through this approach. 

3. The mean convergence for nearby and very distant sources 

Let the origin of the coordinate system be at the observer, who receives a light signal 
at world time r Q . Let symbols like x and y denote the homogeneous Friedmann- Robertson- 
Walker (FRW) coordinate distance in the limit when space curvature has the negligible 
value measured by the microwave background observations (with WMAP being the latest, 
Bennett et al 2003). Suppose the light was emitted by a source at coordinates (x s , y s , 0) and 
passed through en route at world time r a galaxy group with its center at position x (i.e. 
cr = cr — x). If the group is an isothermal sphere, its potential function will have the form 



Further, assume that the middle ray of a pencil beam skirts by the center of the sphere at 
physical distance b (< R), equivalent to coordinate distance y = b/a(r), where a(r) is the 
Hubble expansion parameter at time r. The light is deflected inwards by an angle 



where in Eq. (6) and the rest of this work the speed of light in vacuo is set to unity. Note 
that for b <C R, ip reduces to the familiar constant ip = 2irGM/R = Ana 2 . 

As a result of the lens, the source appears to us to be at position (x s , y' s , 0), where 



The geometry of the situation requires that 




(5) 





y s = y' s - ip(x s - x). 



(8) 
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The apparent size of the source in the radial direction is affected by an amount in accordance 
with dy' s — dy s = (x s — x)dip. The apparent size in the azimuthal direction, however, is larger 
than the true size by the ratio y' s /y s - Hence the effect of the scattering is to increase the 
average angular size 9 of the source by the fractional amount 

e'-e 

" - — < 9 > 

which by virtue of Eqs. (7) and (8) may be written as 



2(l + z)x s \ b db 

where in Eq. (10) and the rest of this paper the expansion parameter at epoch tq is set at 
ao = 1, and z = z(x) is the redshift of the lens. This leads to a fractional increase in the 
energy flux carried by the beam, by the amount 2i]. 



The function z(x) may be obtained by inverting the equation 

1 f z dz' 



x = c(r -r) = — I (11) 



with 



The result is: 



E(z) = = [Q m (l + zf + (1 - Q m )]V2. (12) 

-no 

z(x) = H x + -ft m # V + • • • , (13) 

i.e. an expansion of z(x) as a power series in x. 

In reality the light beam intercepts a random and homogeneous distribution of galaxy 
groups along the way. To work out the total convergence under the scenario of infrequent 
multiple scattering, we first write down the probability of the arriving light having encoun- 
tered a group at position between x and x + dx, and impact parameter 1 between b and b + db, 
as ndV. Here n and dV are respectively the number density of groups and a cylindrical vol- 
ume element, both for the epoch r, i.e. 

ndV = n(T)a(r)dx x 2nbdb = 27m (l + zfdxbdb. (14) 

where n(r) = n(z) has the functional form n(z) = n (l + z) 3 due to Hubble expansion. In 
Eq. (14) the properties of the galaxy groups were assumed not to change with redshift. The 



1 Strictly speaking b should be replaced by the distance of closest approach of the original undeflected ray, 
but in the weak lensing limit where tpx <§C R this distinction is unimportant. 
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validity of this statement is restricted to z being beneath some maximum redshift. Above the 
maximum, the evolution of groups must be taken into account. This limiting value of z seems 
to be at least z ~ 0.5 from direct observation of groups (Jones et al 2002). From dynamical 
time considerations, a group cannot evolve significantly within a timescale < Rja m 3 x 
10 17 s, corresponding to a redshift exceeding z — 1.5. 

We proceed to calculate the expectation value of rj, which is given by 

(77) = ^^i]n5V = J rjndV (15) 

with the integration performed over an entire cylinder wherein lenses are present. By means 
of Eqs. (6), (10), (14), and (15), we performed the integration over b to arrive at 

(77) = 2n 2 GMn H dx [1 + z{x)] ^ ~ ^ , (16) 

JO x s 

where the assumption is for the furthest lens to have a present epoch distance of a^Xf = Xf 
(beyond x = Xf the Universe is too smooth to accomodate galaxy groups as we know them 
today), and that 6 min <C R. In obtaining Eq. (16) use was also made of the definite integral 

1 r R fb\ 

— / arccos — } db — 1 



R Jo \R 

In Eq. (16) we incorporate the contribution to (77) from all types of (galaxy group) 
mass spheres, each having its own mass M and number density n , with ^n M = p c f2 groups , 
where 

* = t§ (17 > 

is the critical density. We can recast the expression for (77) as 

(V) = ^groups^o 2 dX I 1 + <X)] {X " ~ X)X , (18) 

^ JO x s 

which is the mean angular magnification of any small patch of sky randomly located at 
coordinate position (x s ,y s ,0). There are two limiting cases when Eq. (18) simplifies. The 
first is x s = Xf, corresponding to sources embedded within the clumpy environment of the 
near Universe (e.g. Type la supernovae). In this case 

(77) = l -n gTOVLVS Hlx 2 s (\ + ^H x s + iUl m Hlx 2 s + • • (19) 

The 2nd is x s ^> Xf, corresponding to very distant sources which emitted light at a time 
when the Universe was smooth, the cosmic microwave background (CMB). Here we have 

(77) = ^ gr °u P s#o4 + ^HoXf + hl m Hlx) + ■ ■ ^ . (20) 
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Note the absence of any x s dependence in Eq. (20). If corrections due to the finiteness 
of Xf/x s are desired, we remark that the lowest order of such terms equals an additional 
— 2xf/(3x s ) within the last pair of parentheses. 

4. The exact cancellation between the lensing effect of isothermal spheres 
and the Dyer-Roeder demagnification 

A very interesting result which emerges from the analysis thus far concerns the balance 
between beam convergence by the clumps and divergence within the subcritical density 
'voids' between clumps. Take for instance a supernova source, the weak lensing of which is 
described by Eq. (19). If the 'voids' are completely matter-free, i.e. ft g = Vl m in Eq. (19), 
one would attain maximum (77). Yet, in this same limit, the demagnification e of the source 
caused by propagation through the 'voids' will be that of the Dyer-Roeder 'empty beam' 
(Dyer & Roeder 1972), viz. 



and E(z) as defined in Eq. (12). The remarkable fact is that an expansion of e in power 
series of x s with the aid of Eq. (13) gives the same series as Eq. (19) with Vt g = Vt m . 

This development highlights the advantage of an analytical approach. In an earlier 
work, Weinberg (1976) considered an Q\ = Universe where all the matter is clumped into 
point masses, and found that the net magnification is still controlled by the Q m parameter 
alone - as if space remained homogeneous. The present conclusion reinforces Weinberg in 
the more realistic context involving clumps of finite size within a Universe of arbitrary 
and Q m . 

The correspondence between a 100 % clumped Universe and the fully homogeneous 
Universe, as established above, remains in place even if the clumping is not 100 %. This 
was shown in a separate paper (Lieu & Mittaz 2004) where we also adopted a more unifying 
method, using the Sach's optical equations to handle the two opposing effects under one 
formalism. A toy model on the physics of this section is given in Appendix A. The formal 
treatment of average magnification in an inhomogeneous Universe is provided by Kibble & 
Lieu (2005), where the diversity of averages appropriate to different modes of observations 
and data analysis methods will be calculated and discussed. 

5. The standard deviation - convergence fluctuations 

Like the two-point correlation function in galaxy count analysis, the first step towards 



e 




where 




(21) 
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an expression for 5i] is to let rii be the number of groups with their centers lying within a 
volume 5V labelled positionally by an index %. Provided SV is sufficiently small that there 
is no appreciable chance for the centers of two groups to be both inside volume i, then = 
or rii = 1, i.e. 

n\ = n, (22) 

Moreover, because the distribution of groups in space is a Poisson process, and the location 
of each group does not affect that of another, rii has the following properties concerning its 
averages: 

(m) = nSV, (nj) = n5V, (run,) = (n i )(n i ) = (n5V) 2 . (23) 

For each volume there is a corresponding contribution to the convergence. If the i th cell is 
occupied, we will have 

as the fractional change in the angular size of the source due to its light passing by this cell 
(ip' = dtp/db). The total effect is given by a summation along the light path: 

V = ^2 n iVi- (25) 

i 

Taking the average, we have 

(v) = J2( n ^ = J2 n6Vr l*- ( 26 ) 

i i 

i.e. one obtains Eq. (15) for (77). 

Continuing towards the variance, we need the average value of 

= [ Yl HiT]i J = UiT] i + Yl n i n iViVj, (27) 

\ i / i i^j 

which is given, after taking into account Eq. (23), by 

( V 2 ) = J2(n5V) V ^ + J2(^V)\V 3 - (28) 

The variance is now computed in accordance with its definition: 

(s v r = (v 2 ) - (v) 2 = Y> 6V - ( 29 ) 
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where use was made of Eq. (26) and the fact that the (jj) 2 term cancels unless i — j. 
Finally, we note that for small enough 5V, (n5V) 2 <C n5V and may be ignored, so that 
(Srj) 2 = J2inSVri 2 . If the summation over % is cast in integral form, we will have 



rfndV 



87r 3 n cT 4 



In 



R 



dx 



[x s — x)x 

Xe 



(30) 



where in going towards the last expression use was made of Eqs. (10) and (14) with the full 
deflection angle from Eq. (6). 

Note that the standard deviation differs from the mean (rj) in two distinct ways. First, 
while (77) oc ^groups; i^v) 2 1S more complicated, being oc n a 4 . Thus (Si]) 2 is much more 
governed by the specific properties of the type of isothermal spheres in question. Second, 
unlike (77), the integration towards (Sr/) 2 can be performed exactly for all values of Xf and 
x s . When this is done, and the variance for groups of various radii and velocity dispersions 
are summed, the result is 



8tt 



(5 V y = —n xj 1 



2x s 5x 2 J 4^ 



In 



R-i 



7T Z 



(31) 



where pij is the probability of finding a group with velocity dispersion Oi and radius Rj. 



6. Application to Type la Supernovae - a check against numerical simulations 

Utilizing the full statistical properties of the 1,168 galaxy groups observed during the 
ESO survey (Ramella et al 2002, section 1 and Figure 1), the summation in Eq. (31) was 
found to have the value 



In' ^ 



7T^ 



2.90 x 10~ 10 



(32) 



for 2 b m i n = 10 kpc and under the convention of unit speed of light. The equation, together 
with Eqs. (2) and (3), permit an evaluation of the quantities (77) and 5rj as given in Eqs. 
(18) and (30). 

We now apply the results to the astrophysical phenomenon of Type la supernovae 
(SNla), where Xf = x s (see the comment after Eq. (18)), both being < a few Gpc, and 
Eq. (18) reduces to Eq. (19). At redshifts Zf — z s — 1 and 1.5, the values of ((77), 5rj) are 



2 At such a small minimum impact parameter the lensing may have become strong, so that readers can 
legitimately question whether the present calculation, which is based on the assumption of weak lensing, still 
applies. The answer is yes, and is explained in section 8 and Appendix B 
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respectively (0.031, 0.036) and (0.060,0.055). Note in each case 5r) is at most ~ (77) (the 
same applies also to the CMB - see below). This means in the context of our paper there 
is no concern over the probability distribution of 77 becoming asymmetric (skewed). In fact 
the distribution is truncated only at the tail end of the left wing: because from section 4 it 
was proved that (77) exactly equals the demagnification percentage at the voids, a negative 
excursion of 77 from the mean value of (77) by one standard deviation (77) — 5rj does not render 
the source fainter than the rigid lower bound set by the Dyer-Roeder beam. 

In fact, our value of 2Sr] ~ 6 % for the brightness fluctuation at z — 1 compares well 
with that of cosmological N-body simulations by Barber (2000), who obtained 7.8 % (this 
rather large St] may well be the reason why Barber (2000) found an average amplification of 
-3.4 % when our analytical calculation gives zero). We also quote, for completeness, another 
result for 5r/ from an earlier simulation of Wambsganss et al (1997), whose value was 4 %. 
The fluctuation is nonetheless still marginal for the purpose of testing against SNla data, 
wherein the brightness dispersion is at the 0.35 magnitude level (Barris et al 2004, Tonry et 
al 2003). 

7. Application to the CMB - should space be as flat as observed? 

Observations of the CMB are much more accurate than those of SNla, and can be used 
to conduct a useful test of the standard cosmological model. For the CMB at z s m 1000, and 
assuming that the galaxy groups populate the near Universe (without significant evolution) 
to z = Zf we have, from Eq. (18), (77) = 0.099 and 0.157 respectively for Zf — 1 and 1.5. 
Also, from Eqs. (2), (31), and (32), we find that Srj = 0.093 at z f = 1 and 0.133 at z f = 1.5. 

In the present work we focus on 5i], the dispersion in the angular size of CMB hot and 
cold spots. Previous attempts in estimating the extent of this effect involved ray tracing 
codes (e.g. Pfrommer 2003), concerning which it is not entirely clear what galaxy groups 
properties were assumed, and how well they correspond to observations. The quantity Srj 
has the meaning of a standard deviation in the angular size of sources positioned along 
independent sightlines that sample the variation in the spatial location of different sets of 
groups. The separation a between such sightlines is ~ the angular diameter of a typical 
group, with present physical radius ~ 3 Mpc (see Eq.(4)), placed midway (in comoving 
FRW distance scale) between us and z = Zf, i.e. a ~ 12.4 arcmin for Zf = 1, and 9.4 
arcmin for Zf = 1.5. Hence if a spherical harmonic in the CMB TT power-spectrum has 
mean angular size = n/l less than or of order a, the dispersion in 9 will simply be 

59 = 95r] (9 < a, coherent scattering). (33) 

If on the other hand the structure is large, and it contains N ~ 9 2 /a 2 subregions magnified 
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independently by distinct lens combinations, then intuitively 



59 = 



65r] 



(9 > a, incoherent scattering), 



(34) 



although Eq. (34) may also be established firmly by a mathematical proof, which is provided 
in Appendix C. In both cases 59 is to be added in quadrature to the intrinsic dispersion in 
the structure sizes, which are due to density perturbations at decoupling. The dividing 
line between the two cases, a ~ 10 arcmin, is to be compared with the value of 20 arcmin 
advocated by Bartelmann & Schneider 2001 (BS). In fact, if the essential part of our Eq. 



where L — Xf, is converted into the notation of BS, we will make the substitution n Q — * 1/ R? 
(BS assumed L/R lenses per distance L), R — > 1/fc, L — > w, and a 2 /c 2 — > 2$. It will then 
be clear from Eqs. (33) through (35) that 59 is the deflection angle dispersion a{4>) of BS 
for both coherent and incoherent scattering. The key difference, however, is that because for 
the galaxy groups being considered 5rj ~ 0.1, some five times higher than the corresponding 
value in BS, one expects significant broadening of the spherical harmonics. As can be seen 
in Figure 2, such a behavior is inconsistent with the WMAP data. 

An immediate check concerns whether the galaxy group properties we derived from the 
ESO survey are representative of the truth. There are two aspects open to critique. Firstly, 
the cutoff radius given in Eq. (4) is determined from the virial mass M and the dispersion 
velocity a, the former of which (i.e. M) is only an inferred quantity. It is entirely possible 
that in reality we have a mean cutoff radius i?<3 Mpc (i.e. the total mass of a typical group 
is <C its virial mass), but the total amount of matter within galaxies and groups constitute 
an figroup which still satisfies Eq. (3). This would involve a smaller M and R for each group, 
but a larger number density of groups than the value of Eq. (2). In fact, from an earlier 
ESO survey Ramella et al (1999) estimated that if observational selection effects were taken 
into account n group could reach 4 x 10~ 3 Mpc -3 (see Figure 6c of Ramella et al 1999, where 
the number density plotted should be multiplied by h 3 ), i.e. an increase from Eq. (2) by a 
factor of ~ 25. To remain in compliance with Eqs. (3) and (4), then, the cutoff radius R 
(and hence mean angular size a) of the groups must both decrease by the same factor, to 
R ~ 120 kpc. Returning to Eqs. (31) and (32), we see that 5i] becomes somewhat higher, 
but the real difference comes from the \/N reduction factor for incoherent scattering at a 
given CMB spot size, Eq. (34), which is now 25 times more severe. As a result of these 
modifications, the net broadening effect on the primary acoustic peaks is kept drastically 
in check - the theoretical TT power spectrum is no longer distinguishable from that of the 
standard model (the solid line of Figure 2). 



(31), viz. 



(ST,) 2 



n a 



,4 L 3 



(35) 
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The second way of resolving the observational conflict is to question the legitimacy in 
our assumption that all the observed groups are virialized systems with the density profile 
of an isothermal sphere. Strictly speaking, the only relatively secure candidates are those 61 
X-ray emitting groups in the ESO survey (Ramella et al 2002), for reasons already explained 
in section 1. Since these objects constitute « 5 % of the total sample, it could be argued that 
the actual value of 5r] should involve a reduced number density, no — * 0.05no, in which case 
5rj would become « 2 %, and one recovers the minimal distortion advocated by BS, with a 
resulting TT power spectrum that differs negligibly from the solid line of Figure 2. With 
this reduction of 5n, however, the SNla brightness dispersion estimate becomes 0.016, closer 
now to the result of Wambsganss et al (1997) than Barber (2000). One question worthy 
of consideration is, even if most of the groups are not isothermal spheres, they should still 
possess some intra-group gravitational potential, so that the present undertaking of ignoring 
altogether their influence on light may not be justified. Effectively the problem is resolved 
under this scenario by assuming that groups are nothing more than a collection of individual 
galaxies which act as a large number of tiny incoherent lenses, thereby minimizing the size 
dispersion of large light sources. 

8. Generalization of sections 3 - 5 to clumps of arbitrary internal density profiles 

The development of sections 3-5 may readily be extended to the case of mass clumps 
with an arbitrary internal density distribution, with the deflection angle ip{b) of Eq. (6) 
having the general form 

m = 2 f = 4 f" ^Wcosad., (36) 

J-oo r Jo 

where r = y/ x 2 + b 2 = bseca, and the mass m(r) within radius r can be an arbitrary function 
of r. In terms of the density p{r) at r, where 

^ = WpM, (37) 

the fractional weak lensing angular magnification reads 

= 2Gxi(x s - xi) r°° Anrp(r)dr 
V ~ (l + zi)x, J b v^F' 

By integrating over the probability element for randomly located clumps, dP = ndV of Eq. 
(14), and taking account of the fact that 

f°° 4irrp(r)dr , x 

bdb / _^_L = M, (39) 
o Jb Vr 2 - b 2 
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the total mass of each clump, one again arrives at Eq. (18) for the average value of r\. 
Thus the conclusion in section 4 of zero net magnification relative to that in a homogeneous 
Universe of the same mean density is valid for lensing clumps with any profile p(r), and is 
not contingent upon the state of inhomogeneity of the Universe. Moreover, in Appendix B 
we shall prove that, as long as the meaning of the average is appropriately defined, both 
this conclusion and Eq. (18) for (77) remain unchanged even when the light rays are allowed 
to pass through clumps at sufficiently small impact parameters where strong lensing effects 
must also be included. 

In fact, the most general treatment of the problem is to be found in Kibble & Lieu 
(2005), who showed using a more powerful (vierbein) formalism that even under a broader 
range of circumstances than those enumerated above (a) the Universe may still be regarded 
as homogeneous concerning the mean properties of its propagating light, and (b) Eq. (18) is 
the correct expression for (77). One advantage of Kibble & Lieu (2005) lies in its utilization 
of the expansion 9 of a small ray bundle, which is an additive quantity irrespective of the 
strength of the lensing. 

Since the averaging procedure (rj) = J r/dP is applicable to the full range of b, large and 
small, it follows from the derivation of section 5 that the same statement may be made about 
the convergence fluctuation, viz. it too is always given by the expression (Si]) 2 = J rfdP = 
f rfndV . Thus, from Eqs. (38) and (14) we deduce that 



(Srj) 



2 



[ V 2 ndV = 87in G 2 H d Xl (Xj Xl)Xl f°° bdb [ 

J J0 L x s J Jo \_Jb 



4nrp(r)dr 
Vr 2 - b 2 . 



2 



(40) 



for spherical clumps having any p(r). 
9. NFW profiles 

With the availability of Eq. (40) the way is paved for investigating convergence fluctu- 
ation due to the intra-group matter distribution having a profile p(r) which differs from the 
isothermal sphere. In particular, attention is devoted to the increasingly employed Navarro- 
Frenk- White or, NFW, model (Dubinski & Carlberg 1991; Navarro et al 1995, 1996, 1997), 
which involves a function of the form 



*> = drh (41) 



between r = and r = R, where 5 C and r s are respectively known as the overdensity factor 
and scale radius, p c is as given in Eq. (17), and R is the virial radius, related to the r.m.s 
dispersion velocity a by 

R - tto- < 42 > 
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with E{z) as in equation (12) and a as a velocity measured in the frame of the cluster, 
i.e. cr ~ (1 + z). For z < 1.5, (1 + z)/E{z) is a constant. Hence R may also be regarded 
as constant, apart from intrinsic evolution effects which we already argued is without the 
support of evidence in the case of groups. Moreover, the scale radius r s is ~ 26 % of the 
virial radius, or 

— = 3.846 = c (43) 

r s 

(see Carlberg et al 1997), and the overdensity 5 C depends on the ratio c via the equation 

- 200 * (44) 



3 ln(l + c) 



l+c 



Since we are considering the lensing effects of nearby groups, the modification to R by the 
function E(z) is ignored. 



From Eqs. (40) and (41) one arrives at 



(Srj) 



2 



64vr 



R 



3x f 3a; 2 /. 
1 -2i + 55)' (45 > 



\ ^min / \ R ) \ ^mi 

valid for any 6 min except b m \ n <C r s , where the quantity M is given by 

M = 5 cPc r 3 s . (46) 

For galaxy groups with a « 270 km s _1 we obtain, from Eqs. (42) to (44), the following NFW 
model parameters: R = 668 kpc, r s = 174 kpc, and 5 C = 4.83 x 10 4 . Substituting these 
values into Eq. (45), we computed 5r] assuming a group density no such that n^M = Q g , 
where M is the total mass within the virial radius R, and Q g is still given by Eq. (3). 
This involved assigning Hq a value « three times higher than that of Eq. (2). Moreover, by 
setting 3 6 min « 100 kpc, we then find that 5r] is about 1.5 times higher than its corresponding 
value as determined by Eqs. (2), (31) and (32) for isothermal sphere density profiles, when 
the same distances Xf and x s apply to both cases. Thus one can safely declare that the 
conclusions of sections 6 and 7 also hold for NFW density profiles. 

10. The absence of convergence fluctuation in the CMB due to clusters of galaxies 

With all the physical and mathematical prerequisites in place, we are ready for the 
clincher test, which utilizes clusters of galaxies - extensively observed systems with well 



3 We conservatively avoid setting 6 m i n too low, because the performance of NFW profiles in the very inner 
parts of galaxy groups and clusters is controversial. 
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determined properties - as gravitational lenses. Nearby clusters have a number density at 
the present epoch (Bahcall 1988) of 

n « 10" 5 /i 3 Mpc~ 3 (47) 

at h — 0.71, and a mean velocity dispersion in the cluster frame of (Struble & Rood 1991) 
of 

a « 1000(1 + z) km s _1 . (48) 
The latter corresponds, by the low z version of Eq. (42), to a mean virial radius of 

J R^2.12Mpc, (49) 

and hence a mean virial mass in the range 

M = ^^p c R 3 w 1.12 x 1O 15 M . (50) 

Moreover, there is ample evidence that clusters do not evolve significantly put to Zf ~ 1 
(Hashimoto et al 2002, 2004; Maughan et al 2003; Younger, Bahcall, & Bode 2005), so that 
once again for our present purposes n and R may be treated as constants. 

Combining equations (40) through (44) with equations (47) and (49), followed by a 
numerical integration of equation (40) with 6 min = 100 kpc, x s ~ 14.02 Gpc to the CMB, 
and Xf ~ 4.41 Gpc to the farthest lens (set at Zf = 1.5), one obtains 5i] r* 10 %. Given that 
nearby clusters have comparable angular sizes as the galaxy groups (both having R rj 2 - 3 
Mpc), i.e. the effect of incoherent lensing is already calculated in section 7, the implication of 
this result on the TT primary acoustic peaks is again an expected smearing of the gaussians 
as depicted in Fig. 2 (ignoring a slight skewness in the convolution), which is contradicted 
by observations. 

Alternatively one could return to the isothermal sphere model, with p(r) = Mj (47rr 2 i?) 
cutting off at r = R, and the values of R and M as given by equations (49) and (50). Then, 
with the same source-lens distances as in the case of the NFW profile, equation (40) gives 
also 5i] rj 10 % at b min = 100 kpc after a numerical integration, i.e. our conclusion of a 
conflict between prediction and reality remains. 

11 Summary and conclusion 

The all-sky convergence fluctuation of light by galaxy groups, with properties inferred 
from an extensive survey of groups, is computed. When applied to Type la supernovae 
brightness - a situation where the effects are generally hard to measure - the results were 
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found to be in broad agreement with the predictions of earlier authors. When considering the 
CMB, the consequence of a sizable dispersion in the angular diameter of the temperature 
fluctuations is more serious, because it does not correspond to the observational reality 
represented by the latest WMAP data. 

There are at least two ways of overcoming the difficulty, both of which point to the 
possibility that in reality galaxy groups have by far not yet developed to their virial masses 
and virialized potential. The first one advocates the existence of many more groups than the 
directly detected ones upon which the analysis in this paper was based. The total number 
density remains in accordance with the published correction for selection effects (Ramella et 
al 1999), but the mean mass per group is substantially reduced. The outcome is a great deal 
more incoherence in the lensing of the primary acoustic peak structures, and any additional 
dispersion in a spherical harmonic t becomes unobservable. In the second resolution we 
questioned the actual fraction of the groups within the ESO survey sample which have fully 
formed isothermal sphere or NFW profiles. If this is estimated by counting only the X-ray 
emitting groups, the number of such groups is so small that their effect on the CMB will 
again be beneath the margin of detectability. 

Finally, our predicted all-sky variation in the size of the CMB acoustic peaks based upon 
the lensing effect of clusters of galaxies was compared with WMAP data. This raises the much 
more serious question on why convergence fluctuations in the size of the second CMB acoustic 
peak due to clusters are absent. The 'escape route' arguments of the previous paragraph, 
which may have worked in the case of galaxy groups, are no longer so compelling in the 
present context, simply because rich clusters are very well studied systems. Nevertheless, one 
could still contemplate a less profound consequence, viz. perhaps the NFW and isothermal 
sphere profiles are not a good description of clusters after all. Given, however, that the range 
of impact parameters we used to derive the value of Srj is optimal for the performance of the 
profiles, this does not seem to be a sensible way out. We therefore end with the startling 
conclusion that the large scale curvature of space may not entirely be an initial value problem 
related to inflation. The absence of gravitational lensing of the CMB point to the possibility 
that even effects on light caused by wrinkles in the space of the late (nearby) Universe have 
been compensated for, beyond some distance scale, by a mechanism which maintains a flat 
geometry over such scales. 

Authors are indebted to T.W.B. Kibble for his independent re-deriving and cross- 
checking of the mathematical formulae in the paper. 
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Appendix A - Heuristic model to illustrate the coupling between the influence 
of clumps and voids on light propagation 

We provide a simple way of gaining an insight into why the geometry of space as revealed 
by light is, to the lowest order of approximation, independent of the state of inhomogeneity 
of the incipient matter. For the purpose it is only necessary to use non-expanding Euclidean 
space as starting point, i.e. let space be uniformly flat and empty, so that the propagation 
of light is governed by null geodesies in a Minkowski background metric. Clumps may be 
introduced without changing the average properties, by placing at random locations spheres 
of total mass M and internal density profiles p(r) extending to a radius R, with the matter 
for each sphere being drawn evenly from a concentric cavity of radius R\ ^> R. The picture 
is then rather akin to the 'swiss cheese' model. Rays passing at impact parameters b > R\ 
behave as if the relevant clump does not exist. At R < b < R\ a ray is deflected inwards by 
the angle ip — AGM/b because of the clump, and outwards by the angle 



AGM 
^1 = 7— 



(Al) 



because of negative mass distribution in the void. The net inward deflection ip = ip — ipi \s 
then given by 

if, = ^-(R\ - b 2 )l for R < b < R u (A2) 
R{b 

which vanishes at b = R±, thereby satisfying the continuity requirement. Lastly, at b < R 
the clump itself deflects rays inwards according to Eq. (36). With the void included, we 
have 



, AG f w/2 , , , AGM 

ip = —— / m(r)coso;ctQ; 

o Jo b 



(A3) 



where r = feseco; and m(r) is the mass within radius r (m(R) = M). Again, the two values 
of if) at the b = R boundary match. 

Next, let the source be at an infinite distance away, and the clump- void system be 
centered at distance x from the observer. The angular magnification formula of Eq. (10) 
reduces, in the present context, to 

Note that when light transits the void region (i.e. R < b < Ri) there is demagnification - this 
is the Dyer-Roeder effect. The average weak lensing magnification is obtained by integrating 
i] w.r.t. the probability element appropriate to a random distribution of clump- void systems 
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of density no, viz. dP = lim^hdhdx . If the clumpy region spans the range x = and x = L, 



(rj) = 47fGn / xdx 
Jo 



bdb 



47rrp(r)dr 3M 



Vr 



b 2 



R\ 



1-^2 I ^ 



0, (A5) 



where Eq. (39) was used to calculate the first term in the square parentheses. Thus, as long 
as the light signals pass through a sufficient number of clumps and voids, the average source 
size does not deviate from the Euclidean benchmark. 



Appendix B - on how to average the magnification if strong lensing is included 

To establish the correct averaging procedure it is only necessary to consider one fore- 
ground clump, placed at the center of the field of a large, circular, and uniformly illuminated 
background source. Take a small annulus of emission, the undeflected and actual light paths 
connecting it and the observer pass by the clump at impact parameter & an d b respectively. 
In the weak lensing limit, it is unnecessary to distinguish b from b . In the absence of the 
lens, the solid angle subtended by the annulus at the observer is 

,~ /-, , 2 2nbdb 

dw = (l + *j) 2 — . 

x i 

In the presence of the lens dw is increased by the amount 47i(l + zi) 2 r]bdb/x 2 where, according 
to section 8, 77 is given by 

with L = xi(x s - xi)/[(l + zi)x s \. 

If the unlensed source occupies (in projection on the lensing plane) a circular area of 
physical radius Bq as measured at z — zi, the average percentage magnification of its observed 
solid angle may be expressed as 



5Cj 1 o nr r^rp(r)dr SGML 



where in the second last equality used was made of Eq. (39), and in the last equality the 
average of i] is defined as 

= slpm 

f Q B 2nbdb 

and matches Eqs. (14) and (15) of the main text with 1/ J 2nbdb replacing n(z)Sxi/(l + z{), 
since in this Appendix the number of clumps within a comoving slice at distances x\ — > x t +Sxi 
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is exactly one. Note that from Eq. (B2) the percentage magnification depends only on the 
mass of the enclosed clump - it is not affected by the details of the clump's internal density 
profile. This is reasonable, because all forms of lensing conserve surface brightness (i.e. no 
overlap of emission regions in the image possible) the extra solid angle the magnified image 
claims is simply given by the amount of outward deflection of the 'boundary rays' of the 
source - a process which relates only to the clump mass - because such rays are too far away 
from the clump to be mindful of its inner structure. 

On the other hand, we know that the above consideration is limited to the regime 
of weak lensing, so the question is whether inclusion of strong lensing modification of the 
integrand in Eq. (Bl), which is inevitable as one reaches the bottom end of the integration, 
would lead to new terms carrying extra parameter dependence to jeopardize our hitherto 
consistent result. 

To see why the answer is no, we must now extend the treatment to accomodate the 
possibility of strong lensing, i.e. the relationship b m b is no longer rigorous enough. Rather, 
it has to be expressed more precisely as 

b = b- 14(b). (BA) 

We may define the reciprocal magnification of the annulus area by the formula 

J=- = ^ (B5) 
\i bdb 

We see that while the average direct area magnification ji 

<"> = ^ 

with asymptotic boundary parameters related by 

, x AGML , x 

B = B - L^j(B) = B — , (B7) 

is divergent because at sufficient low b one encounters the caustic J = which is clearly not 
a weak lensing phenomenon, the average of J = is free from infinity problems. Thus our 
undertaking to evaluate the average reciprocal magnification is not motivated by the need 
to obtain the 'desired' outcome, but to obtain a meaningful outcome. Explicitly 



1 \ J B 2nJbdb B 2 



2 



<^y^^> (B8) 
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where in reaching the last expression use was made of the fact that J is the Jacobian of 
transformation from bdb to b db . Note also from Eqs. (B3), (B6), and (B8) that the 
calculation of all the averages are consistent. 

The average percentage magnification of the image area, including the effect of strong 
lensing of the central rays, is now equal to (J) -1 — 1 = (B 2 — Bq)/B 2 . By Eq. (B7), 
this is just 8GML/B 2 , in agreement with Eq. (B2). We succeeded in proving that the 
formula for (rj) in Eq. (B3), which as explained in the material immediately following this 
equation reflects exactly the same averaging procedure adopted in the main text of the entire 
paper, is appropriate to the full regime of weak and strong lensing. Therefore, the statistical 
cancellation between lensing magnification and the Dyer-Roeder beam, sections 4 an d 8, is 
a robust conclusion not contingent upon the strength of the lensing. Our Eq. (40) for the 
variance (Sr/) 2 is, for this reason, also valid at arbitrarily low impact parameters. 

Those interested in 'seeing beyond the mathematics' to actually understand how the 
numerator integral for (J) in Eq. (B8) yields -kB 2 despite the presence of multiple images 
in strong lensing should consider in detail how the integration is carried out in b space 
Specifically, in Eq. (B4) b does not change monotonically with b. As the latter decreases 
there comes a place, say b = b±, at which the corresponding 6 = ^10 = 0. After this, initially 
for b < bi, bo turns negative. For any non-singular central density function p(r), however, 
ip(b) — > as b — > 0, so when b decreases further from b — b\ towards b = a point b = b 2 
will come at which 6 reaches a minimum of b = 6 20 , thereafter increasing again towards 
b = 0. Obviously J = at both b = bi and b = b 2 - Within the strong lensing regime, then, 
there are in general three images of any given source pixel, as is well known. The first one 
corresponds to b > b±, and is on the same side of the center (on-axis position) as the source 
pixel. The second and third images correspond to b 2 < b < bi and < b < b 2 respectively, 
and lie on the opposite side. For the first and third image, J > 0; for the second, J < 0. 
The J-integral may therefore be separated into three, viz. 

pB pB pbi pb 2 

/ 2nJbdb= / 2jrJbdb+ / 2nJbdb+ / 2nJbdb. (B9) 

JO Jbx Jb 2 JO 

Changing variable to y = \ b 1 and applying Eq. (B5), one may see how the integration takes 
place in b space, with the sign of J for the three images assuming particular importance: 

pB pB /■ 1 620 1 p\b2o\ 

/ 2nJbdb = / 2nydy - / 2ixydy + / 2%ydy = irB 2 , (510) 
Jo Jo Jo Jo 

in agreement with Eq. (B8). 

Appendix C: On the convergence fluctuation of a large source 
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Like the problem in Appendix B, the necessary theorem here may be established by 
considering only one slice of redshift. Within this slice let us introduce transverse coordinates 
y and suppose that for one clump at yi the value of 77 for a small bundle of light rays at y is 
given by r)(y) = f(y — yi)- If the number of clumps per unit area is n, the average of rj will 
be 

m=*J f(y-yi)d%. (Cl) 

If the averaging is performed over a sufficiently large area fj(y) = (rj) should become inde- 
pendent of y. Moreover, the equation 

(5 V f = nJ f 2 (y)d 2 y. (C2) 

gives, by section 5, the variance in (77). 

It's also useful to define the correlation function £(y) via 

(5 v n(y) = n J f(y')f(y + yW- (C3) 

Clearly, £(0) = 1, and falls off with y on a scale characteristic of the size of a clump. 
We can therefore write 

Ai = J md 2 y- (C4) 

as definition of the area A\ of a clump. 

To calculate the mean and variance in the magnification of a background source of 
projected area A at the lensing plane, let us first consider the effect of all the clumps in 
a really large area S ^> A. The probability that within S there are exactly M clumps 
at the positions y±, . . . , y^, say, within small areas d 2 y\, . . . , d 2 yjy, is given by the Poisson 
distribution 

d 2N P = e- s n N \[d 2 y r (C5) 
The total increase in the area of A due to lensing by these clumps is 2r/ A A where 

Va a= / d 2 yJ2f(y-yj)- (C6) 

3=1 

To find the average, we integrate (C6) over the measure (C5), and sum over Af. The result 
is 

( VA )A = / d 2 ye- s Tn / ^ ■ ■ ■ ^ E 



- 21 - 



Note that the J\f\ factor is needed to compensate for multiple counting of permutations of the 
clumps. Also, as already explained in Appendix B, this averaging procedure has correctly 
taken into account strong lensing effects as well. Thus we obtain 

(V), (C8) 



1 r 



where in reaching Eq. (C8) we observe that each J s d 2 Xi term in the summation over j in 
Eq. (C7) gives the same contribution: one of the integrals over Xi yields, by Eq. (CI), (rj) /n, 
and the others altogether yield S^ -1 . 

Concerning the average of rj\, we now have to take the square of Eq. (C7) and integrate 
over the measure of Eq. (C6), i.e. 

( Va )a* = / #y / d 2 g e -» s e Tn / . . . d 2 m E f(y- y^tf - ^ 

Ja Ja tf=l A,lJ S j,k=l 

In Eq. (C9) it is necessary to separate the terms with j = k from the others. All terms with 
with j 7^ k are equal, totalling 



00 n-^- 2 r r 



i.e. these terms add up to the very quantity that needs to be subtracted from (^)A 2 to 
obtain the variance. The variance is then simply given by the sum total of all the j — k 
terms. We arrive at 

(5 VA ) 2 A 2 = e~» s E Jj^JyS^ 1 J/y //¥ J s d% n f(y- ft)ftf - ft). (Cll) 
By Eq. (C3) we have, finally, 

(5 VA fA 2 = (5 V f [ d?y [ d^i{y - V). (C12) 

JA 

Defining N = A/Ax as the ratio of the source area to the clump area, then, so long as 
N is large (i.e. away from the N < 1 regime) y will not be too near the boundary, and the 
integral over y 1 will give Ai, while the remaining integral over y is = A. Thus, apart from 
small (boundary-effect) corrections, 

(SvAf = (Srjf^ = {CIS) 
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Note added in Proof (dated 14th June, 2005): We are now aware of the paper by Scranton 
et al (ApJ in press, astro-ph/0504510), which reported the first secure detection of magnifi- 
cation dispersion, or convergence fluctuation, among a large sample of distant quasars. This 
provides the necessary observational evidence that the effect discussed in our present work 
is real, and further elevates the profile of the question on why the CMB exhibits no sign of 
having been lensed. 
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Fig. 1. — The probability distribution of galaxy group velocity dispersion (top panel) and 
the distribution of group mass against velocity dispersion (bottom panel), both derived from 
the ESO survey (Ramella et al 2002). These data have been used to evaluate the quantity 

T,i,jPij a iM R j/ b min) Of Eq. (31). 
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Fig. 2.— The standard FRW model with h = 0.71, Vt m = 0.27, fi A = 0.73 as applied to the 
TT power spectrum measured by WMAP is plotted as a solid line. If in this model 50 % 
of the matter within z ~ 1 is clumped into isothermal spheres with properties given by the 
ESO survey of galaxy groups, the lensing induced convergence fluctuation (5rj ~ 10 %) will 
cause a smearing of the spherical harmonics at and around the primary acoustic peaks by 
an amount in accordance with Eqs. (33) and (34). The resulting theoretical power spectrum 
is plotted as a dashed line - it is an unacceptable fit to the data, with Xrcd = 

2.42 for 38 

degrees of freedom (null probability hypothesis < 0.0001). 



